Photocurrents in nanotube junctions 
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Photocurrents in nanotube p-n junctions are calculated using a non-equilibrium Green function 
quantum transport formalism. The short-circuit photocurrent displays band-to-band transitions 
and photon-assisted tunneling, and has multiple sharp peaks in the infrared, visible and ultraviolet 
ranges. The operation of such devices in the nanoscale regime leads to unusual size effects, where 
the photocurrent scales linearly and oscillates with device length. The oscillations can be related 
to the density of states in the valence band, a factor that also determines the relative magnitude of 
the photoresponse for different bands. 



Since their discovery 0, Carbon Nanotubes (NTs) have 
been the subject of intensive research due to their in- 
triguing electronic and structural properties, and have 
demonstrated great promise for future nano-electronic 
devices 0- However, their potential for opto-electronic 
applications has received much less attention, despite the 
ideal properties that NTs present. Indeed, a desirable 
property of opto-electronic materials is a direct band- 
gap, since it allows optical transitions to proceed without 
the intervention of phonons. NTs are unique in this as- 
pect, since all the bands in all semiconducting NTs have 
a direct gap. Furthermore, the low dimensionality of NTs 
leads to a diverging density of states at the band edge and 
a high surface-to- volume ratio, reducing the sensitivity to 
temperature variations and allowing efficient use of the 
material. Finally, non-radiative transitions can signifi- 
cantly reduce the performance of conventional materials; 
NTs are believed to have low defect density, reducing 
non-radiative transitions. 

These unique properties of NTs are only starting to 
be explored for opto-electronics. Recent experimental 
work|jj has shown that NT field-effect transistors can 
emit polarized light, while illumination of these and 
other NT devices generates significant photocurrent 0,0- 
Such observations have also been made in semicon- 
ductor nanowires0, Q, El- Often, the observed opto- 
electronic effects are due to the presence of a p-n junc- 
tion, with various physical realizations (electrostatically 
defined^, Q, crossed- wire geometry 0, or modulated 
chemical doping §.) Modulated chemical doping of in- 
dividual semiconducting single-wall NTs to create "on- 
tube" p-nj unctions has recently been reported in the 
literature^- These NT p-n junctions serve as excellent 
testbeds for understanding opto-electronics at reduced 
dimensionality. 

Here we show that in simple NT p-n junctions, the 
photocurrent shows unusual features. Unlike traditional 
devices, the photoresponse in the NT junctions consists 
of multiple sharp peaks, spanning the infrared, visible 
and ultraviolet ranges. Furthermore, at nanoscale di- 
mensions, the NT junctions show size effects, where the 
photocurrent scales and oscillates with device length. 

We consider a NT p-n junction under illumination, as 



FIG. 1: Sketch of a portion of the scattering region for the 
(17,0) nanotube p-n junction. The scattering region is at- 
tached to two semi- infinite leads (not shown in figure). The 
incoming light of energy hw is polarized along the length of 
the nanotube, with electric field E and magnetic field B. 



shown in Fig. ^ We use a (17,0) single-wall, semicon- 
ducting NT of radius 0.66 nm, and model the doping as 
in Ref. [lfj. Our tight-binding Hamiltonian with one 7r 
orbital per carbon atom and nearest-neighbor matrix ele- 
ment of 7 = 2.5 eV gives a direct band-gap for this NT of 
0.55 eV. The incident light is assumed to be monochro- 
matic, with polarization parallel to the NT axis. 

To calculate the photocurrent, we use a non- 
equilibrium Green function formalism [llj . in a real- 
space representation. The device consists of an illu- 
minated scattering region connected to two dark semi- 
infinite leads, which are also doped (17,0) NTs. The 
structure of the (17,0) NT corresponds to parallel rings 
of 17 atoms, which are spaced by 0.07 nm and 0.14 nm. 
Each ring forms a layer in our system, with the rings in 
the scattering region labelled from 1 to N. 

Since only layers 1 and N couple with the leads, the 
steady-state, time-averaged current can be written as 
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(This expression assumes that the hole current is equal 
and opposite to the electronic current and that the cou- 
pling of layers 1 and N to the leads is 7. The appropri- 
ateness of this equation for time-dependent Hamiltonians 
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is discussed in Ref. 01) Here, G^ +1 N and Gf Q cross 
the scattering region boundary. To relate these terms 
to G < calculated entirely in the scattering region, we use 
the Dyson equations 13] 7 G< = G^E^ + E^G^ and 

^Gn+i,n = Gn,n^n!n + ^n°nGn,n- (We also note that 
in the presence of electron-photon interactions and the 
potential step in the v-n i unction fl^|. G^ +1 N ^ Gf Q , 
an equality that would normally hold for a system with- 
out inelastic scattering.) 

In the scattering region, G < is determined from the 
equations 



G< = G fl E<G flt 



(2) 



and 



approach is that a model for T,^^ does not have to be 
formulated. (Including electron-hole recombination and 
other multi-photon processes requires going beyond the 
perturbation scheme introduced here, and is beyond the 
scope of this paper. At low light intensities, our calcula- 
tions should provide a reasonable approximation.) 

To derive the function E < (p' 1 ) due to the electron- 
photon interaction, we use the interaction Hamiltonian 



Hel-ph — — A • p 
m 



(6) 



where A is the time-dependent electromagnetic vector 
potential, p is the electronic momentum operator and 
m is the electron mass. Following the procedure of Ref. 
Ha. we obtain 



G R = [EI-H-Z R ]-\ 



(3) 



where H is the Hamiltonian for the NT p-n junction un- 
der illumination. The functions S fl and E < represent the 
interaction of the scattering region with the semi-infinite 
(17,0) NT leads and the incident light. In the presence of 
the electron-photon interaction, these functions depend 
on G < , and the equations are coupled and non-linear. 
To simplify the calculation, we perform an expansion to 

first order in the light intensity: G< = G <0 + G <( - ph \ 
g r = G R0 + G R( P h)^ s < = E <o + E <( P h) and £ jt = 

£_r° _|_ ^.RO/i) ^ wnere the order zero functions represent 

the dark p-n junction. In the absence of an applied bias, 

the only current in the device is the photocurrent 



j( P h) = f£2 i ilE R(1 
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where 



+G™£ <0 G R0 Z R ^G 



Q<(ph) _ QR0ji<(ph) qRO^ I QRO^Riph) qR0y, <0 G R0 ^ 

/• (5) 

Here G R0 = [EI - H - X R0 ] ~\ with the Hamiltonian 
Hq and the self-energy function £™ describing the p- 
n junction without incident light, including the electro- 
static potential. (The Hamiltonian matrix elements are 



tj21,21-1 Tj2Z-l,2i 
Ad 



2 7 cos(^),^ 



H, 



21+1,21 




7, and H ' = —eVi, where Vi is the electrostatic poten- 
tial on layer I). The self-energies E fl0 due to the semi- 
infinite leads are calculated using a standard iterative 
approach [li). and £ <0 = -2fImY> R0 where / is the 
Fermi function. The electrostatic potential and charge 
distribution are determined using a self-consistent proce- 
dure as in Ref. 0. The last two terms in Eq. (jSJ are 
proportional to f(E), while, as discussed below, the first 
term is proportional to f(E — tvjj). As long as the Fermi 
level is in the bandgap and not too close to the conduc- 
tion band edge, the ratio f(E)/ f(E~hoj) is small, so we 
neglect the last two terms in Eq. J5J . A benefit of this 



^<(ph) 
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Pi p P qm G<»(E-Hu>), (7) 



with 



Pirn = < 



— 1 I — m = — 1 , I even 
+ 1 I — m = 1, I odd 

— cos Cjjf) I —m= —1, I odd 
+ cos (-jj) I — m = 1, I even 
otherwise 



(8) 



In these equations, a = 0.07 nm is the smallest ring 
separation, fku is the photon energy, c is the speed of 
light and e is the permittivity of free space. J is the 
circumferential angular momentum, labeling each of the 
AI = 17 bands of the NT. Conservation of J is en- 
forced, as required for light polarization parallel to the 
tube axis^(|. The photon flux F represents the strength 
of the incoming radiation, in units of photons per unit 
time per unit area. Since the photocurrents calculated 
in this work are proportional to F, we focus on the 
photoresponse I/eF. (This definition of the photore- 
sponse might seem awkward at first, having units of 
area/photon. However, the unusual size effects in the 
nanotube devices are more transparent with this defini- 
tion.) Because G <0 = -2/ImG™, T^ h) (E) is propor- 
tional to / [E — hoj). 

Figure [2K shows the calculated self-consistent band di- 
agram for the NT p-n junction for a doping of 5 x 10~ 4 
electrons/C atom. The figure shows the conduction and 
valence band edges of the band with the smallest band 
gap, J = 6, shifted by the local electrostatic potential. 
The band bending is characterized by a step at the junc- 
tion and flat bands away from the junction region 17]. 

The calculated photoresponse for this device at zero 
bias when 128 rings are illuminated is shown in Fig. [5Jd. 
The photoresponse due to bands J = 6,5, 7, 4, 3, 8 (in- 
creasing band gap) is separately plotted in the figure 
(bands 11,12,10,13,14,9 are equivalent). The general 
trend is for the photoresponse for the different bands J 
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FIG. 2: Panel (a) shows the calculated self-consistent band 
diagram for the NT p-n junction. Solid lines are the valence 
and conduction band edges for the J = 6 band. Dotted line 
is the Fermi energ level. Path 1 represents a band-to-band 
transition while paths 2 and 3 represent photon-assisted tun- 
neling. Panel (b) shows the calculated photoresponse for this 
device, plotted seperately for each band J (labelled with the 
value of J in the figure. 



to peak at higher photon energies as the band gap in- 
creases. Because the scattering cross-section decreases 



with fko 



one would ex- 



the (fku) factor in Eq. Q 

pect the maximum photoresponse attained for each band 
to decrease with band gap. Surprisingly, the height of the 
peaks in Fig. Et> does not follow this behavior; in partic- 
ular the response for band J = 8 is actually larger than 
the response for band J = 6. This is a result of the differ- 
ent bands having different effective masses. Indeed, the 
effective mass for J = 8 is about 36 times larger than for 
J = 6, leading to a much larger density of states near 
the edge of the valence band. (The role of the density of 
states will be discussed further below.) 

Path 1 in Fig.[2Ji shows a band-to-band transition with 
the absorption of a photon of energy h\o = 0.6 eV. An 
electron coming from the p-type side of the device in 
the valence band absorbs a photon and is excited to the 
conduction band, where it then continues to the n-type 



terminal. Such a transition is allowed when the photon 
energy exceeds the band gap E g (0.55 eV for band 6). 
Band-to-band photocurrents in the NT device in the left 
lead due to electrons coming from the n-type terminal 
vanish unless the photon energy is larger than the band 
gap plus the potential step across the junction, huj > 1 
eV in Fig. |2Jd. This asymmetry in the currents to the 
left and right terminals leads to the net photocurrents in 
Fig. 03. 

While these band-to-band transitions explain part of 
the photoresponse, a significant response exists at en- 
ergies below the band gap. Such contributions can be 
attributed to photon-assisted tunneling. Paths 2 and 3 
in Fig. [2K show two possible paths for photon-assisted 
tunneling. For a given band J, this process can only oc- 
cur when hio > A,/, where Aj is the difference between 
the asymptotic conduction band edge on the p-type side 
and the asymptotic valence band edge on the n-type side 
(equal to 0.06 eV for the J = 6 band shown in Fig. EK) - 
The photon-assisted tunneling thus turns on at hio = A,/. 
As the photon energy increases above Aj, more states in 
the band gap become available for transport, and the 
photoresponse increases. For the bands with larger band 
gaps, the tail due to photon-assisted tunneling is less im- 
portant relative to the band-to-band peak, since tunnel- 
ing probabilities decrease with band gap. 

The photoresponse of the different bands leads to mul- 
tiple sharp peaks in three different regions of the elec- 
tromagnetic spectrum: infrared, visible and ultra-violet. 
This unusual behavior arises because all the bands in 
the NT have a direct band-gap, which leads to a re- 
sponse over a wide energy spectrum. The separation of 
this wide response into peaks grouped in different regions 
of the electromagnetic spectrum is due to the particular 
electronic band structure of the NT, which has groups of 
bands separated by relatively large energies. The conduc- 
tion band edges for J = 6 and J = 5 (infrared response) 
are separated from those of bands 4 and 7 (visible re- 
sponse) by about 0.6 eV, which arc in turn separated 
by about 0.5 eV from the J = 3 and J = 8 conduction 
band edges (ultraviolet response). Since the photocur- 
rent consists of an excitation from the valence band to 
the conduction band, the separation in Fig. |2b is twice 
these values. 

In conventional bulk junctions, the photoresponse de- 
pends only on the dimensions of the device perpendicular 
to current flow. The nanotube device however shows a 
dependence with length, as shown in Fig. |3| Panel (a) 
in the figure shows the photoresponse of band J = 6 for 
different lengths of the illuminated region. As the illu- 
minated region becomes longer, the peak in the photore- 
sponse increases and also moves to lower energies. This 
behavior leads to linear scaling of the photoresponse with 
the length of the illuminated region, as the integrated 
photoresponse in the inset of Fig. shows. The pho- 
toresponse for three different photon energies is shown in 
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FIG. 3: (a) Dependence of the photoresponse for band J = 6 
on the length of the illuminated region. The inset shows the 
integrated photoresponse, and its linear scaling with length. 
The response for different photon energies is shown in (b). 



Fig. 03. Clearly for hw = 0.4 eV (photon-assisted tun- 
neling) the photoresponse saturates with length, due to 
the fact that the wavefunctions in the band gap decay 
exponentially away from the junction. The response for 
hw = 0.612 eV and hw = 0.7 eV shows a completely 
different behavior, oscillating around a general linear in- 
crease. 

These surprising results can be understood from a sim- 
plified description of the photocurrent in the nanotube. 
To this end, we note that the photocurrent density flow- 
ing to the right lead (i^ ^ = ^| J i R (E)dE^ can be ex- 
pressed as 

i R (E) = -T R Im ( G ™ S <(^) G B0t) (9 ) 

where T R — — 2ilm£^jy. (Eq. © can be obtained 
from Eq. Q by using Dyson's equation r yG^ r+1 N — 
G^ N Z%° N + Z%° N G^ N , the relation G< = -Gf* and 
the condition that the net current must vanish for the 



dark junction at zero bias.) At zero bias, G <0 is purely 
imaginary, and therefore so is T, < ^ ph \ The above equa- 
tion becomes 

i R (E) = -F R £ Im Z< PH) GZ Re G™ + Im G™ Im G™ } . 

k.p 

(io) 

For a given outgoing electron energy E and photon en- 
ergy hw, we have found numerically that the argument 
of the summation is peaked around the diagonal k = p. 
Taking only these contributions, we obtain 

i R W~-r R ^TmTi<™\G%},\'. (11) 

k 

The photocurrent in the NT can thus be understood in 
terms of the excitation of electrons in each layer k along 
the NT (the Im K^ p ' 1 ' term in the above equation), and 
the subsequent transmission to the lead by G™. . 

To explain the linear scaling, we note that for hw > 
Eg and E > E~°° where E~°° is the asymptotic value 
of the conduction band edge on the p-type side of the 
device, there is a section of the NT where band-to-band 
transitions are allowed, and the sum in the last equation 
is dominated by this section of the NT. As the length of 
the illumination region is increased, a longer section of 
the NT is available for band-to-band transitions, leading 
to the linear scaling of the current with length. 

While this explains the linear scaling, there remains 
the questions of the oscillations in the photoresponse 
and the dependence of the photoresponse on the effec- 
tive mass. To address this, we note that lmY,^ ph \E) 
contains terms like Im G^2{E — hw), which are related to 
the density of states at layer k, Dk, through Im G^(E — 
hw) = 7tf(E - hw)D k (E - hw), with f{E) the Fermi 
function. Therefore, the photocurrent is sensitive to the 
density of states at energy E — hw. This explains the 
origin of the dependence on the effective mass, and the 
relative height of the peaks in Fig. [21 

Figure 0| shows the valence band density of states for 
J = 6, calculated at layer 11 for systems with illuminated 
lengths of 24.78, 26.88, and 28.98 nm. (This is the density 
of states near the edge of the illuminated region, which is 
moving further away from the junction as the length in- 
creases). The density of states contains many peaks, and 
as the system size changes, the peaks move in energy. 
Because the propagator | G™ | is sharply peaked at en- 
ergy E~°°, the photoresponse is particularly sensitive to 
the density of states at E~°° — hw. At that energy, the 
density of states oscillates as a function of the distance 
from the p-n junction, as illustrated in the inset of Fig.0] 
This leads to the oscillations in the photoresponse as a 
function of illumination length shown in Fig. [3Jd. The 
oscillation wavelength of the density of states increases 
for energies closer to the band edge, causing the different 
oscillation wavelengths for Hw = 0.612 eV and hw = 0.7 
eV in Fig.EJj. 
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FIG. 4: Density of states at layer 11 as a function of energy. 
Solid, dashed, and dotted lines are for illumination lengths of 
24.78, 26.88, 28.98 nm respectively. The top (bottom) inset 
shows the density of states on the even rings at energy -0.1 
eV (-0.2 eV). 



Although the device presented here is fairly simple, it 
already shows the richness of new phenomena that arises 
in nanoscale opto-electronics. One may envision several 
uses for this or related devices (photodetection, power 
generation, optical communication, etc.), but the impor- 
tant point is that the behavior is much different from 
traditional devices and must be taken into account in 
future device development. 
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